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O ■ ABSTRACT 

o 

f — ■ We make use of interferometric CO and Hi observations, and optical integral-field spec- 

^^ I troscopy from the ATLAS''^ survey to probe the origin of the molecular and ionised inter- 

stellar medium (ISM) in local early-type galaxies. We find that 36±5% of our sample of fast 
rotating early-type galaxies have their ionised gas kinematically misaligned with respect to 
the stars, setting a strong lower limit on the importance of externally acquired gas (e.g. from 
mergers and cold accretion). Slow rotators have a flat distribution of misalignments, indicating 
that the dominant source of gas is external. The molecular, ionised and atomic gas in all the 
detected galaxies are always kinematically aligned, even when they are misaligned from the 
stars, suggesting that all these three phases of the interstellar medium share a common ori- 
gin. In addition, we find that the origin of the cold and warm gas in fast-rotating early-type 
galaxies is strongly affected by environment, despite the molecular gas detection rate and mass 
fractions being fairly independent of group/cluster membership. Galaxies in dense groups and 
the Virgo cluster nearly always have their molecular gas kinematically aligned with the stellar 
kinematics, consistent with a purely internal origin (presumably stellar mass loss). In the field, 
however, kinematic misalignments between the stellar and gaseous components indicate that 
at least 42±5% of local fast-rotating early-type galaxies have their gas supplied from exter- 
nal sources. When one also considers evidence of accretion present in the galaxies' atomic 
gas distributions, ^46% of fast-rotating field ETGs are likely to have acquired a detectable 
amount of ISM from accretion and mergers. We discuss several scenarios which could explain 
the environmental dichotomy, including preprocessing in galaxy groups/cluster outskirts and 
the morphological transformation of spiral galaxies, but we find it difficult to simultaneously 
explain the kinematic misalignment difference and the constant detection rate. Furthermore, 
our results suggest that galaxy mass may be an important independent factor associated with 
© 2009 RAS the origin of the gas, with the most massive fast-rotating galaxies in our sample (Ma' ^ — 24 

mag; stellar mass of wSx 10^" M0) always having kinematically aligned gas. This mass de- 
pendence appears to be independent of environment, suggesting it is caused by a separate 
physical mechanism. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: ISM - 
TSM' molecules — TSM' evolution — stars' mass-loss 
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1 Introduction 



I Dunlap Fellow 



Lenticular and elliptical galaxies, collectively known as early- 
type galaxies (ETGs), have long been thought to be the endpoint 
of galaxy evolution. These systems have surprisingly uniform red 
optical colours and are located i n a tight red sequen ce in an optical 
colour-magnitude diagram (e.g. ISaldrv et alj|2004f) . The origin of 
the colour bimodality that separates the ETGs in the red sequence 
from star- forming galaxies, located in a 'blue cloud', has become a 
central question in extragalactic astrophysics. In order for galaxies 
to fade and join this tight red sequence rapidly enough, it is thought 
that the fuel for star formation must be con sumed, destroyed or re- 
moved on a reasonably short timescale fe.g. lFaber et alj2007h . 

In the paradigm of cold dark matter with a non-zero cosmolog- 
ical constant (ACDM), it is thought that star-formation is shut-off 
in these galaxies by majo r mergers that can funnel the cold gas into 
the centre of the remnant ( lBarnesl2002h . where it may be destroyed 
by an active galactic nucleus (AGN) or be used up quickly in a star- 
bu rst. In massive galaxies, v irial shocks (such as those discussed 
by ISimboim & Dekell 120031) can also create a halo mass thresh- 
old, above which incoming gas will always shock to the virial tem- 
perature, stopping the accretion of further cold gas. Energy input 
from an AGN or starburst may also remove cold gas without a ma- 
jor merger, and keep halo gas from cooling. Galaxies in clusters 
and groups may also have their interstellar medium (ISM) removed 
by harassment, or stripped and ablated away by a hot intra-cluster 
medium (ICM). 

Evolutionary processes like those listed above leave galax- 
ies on the red sequence with little or no cold ISM, and thus 
no star formation. Evidence is mounting, however, that a reason- 
able fraction of ETGs do have cold gas re servoirs and r esidua l 
star-formation. For e xam ple the surveys of Grossi et al.l OOOgh . 
iMorganti eT^ l l2006h and lOosterloo et al.l bOld) find between 44 
and 6 6% of field ETGs have detectable neutra l hydrogen reservoirs, 
while lColbert. Mulchaev & Zabludofll l l200lh have shown that over 
75% o f ear ly-types c ontain dust features in optical images. I Yi et alj 
( I2OO5I) and lKaviraret al. ( 2007) have also shown that at least 30% 
of ETGs have an excess of UV emission, attributable to ongoing 
residual star formation (over and above that predicted for an old 
stellar popula tion from the UV- uptum phenomenon; see lYiluOOa 
for a review). ISarzi et al.l ( 120061) have shown that at least 10% of 
local ETGs have strong Balmer line emission and line ratios indica- 
tive of ongoing star-formation. These studies all imply some ETGs 
must have a cold molecular gas reservoir to fuel this ongoing star- 
formation. 

The molecular gas content of early-types had until re- 
cently been investigated main l y usin g infrared-selected surveys. 
For example iKnapp & Rupenl ( 119961) report a molecular gas de- 
tection rate of 80% for ETGs brighter than 1 Jy at 100 /xm. 
ICombes, Young & Bureaul ( 120070 performed one of the first stud- 
ies of ETGs which was not infrared-selected. They observed 
the SAURON sample of galaxies (a representative sample of 
48 ETGs which had been observed with the SA URON inte- 
gral field spectrograph, see Ide Zeeuw et alj 120021 for full de- 
tails) and reported a molecular gas detection rate of 28% (with 
a CO(l-O) se nsitivity of 3mK T* in a 30 km s^^ channel). 
The survey of IWelch. Sage & Yound (I2OI0I) (which combined and 
added to the results of IWelch & Sagdl2003l : isage & Welclj|200^ : 
ISage. Welch & Young|2007l) was the first published volume-limited 
sample, which found a detection rate of 26% for ETGs within 20 
Mpc (with a similar sensitivity). 

The ATLAS^'' project is a complete, volume-limited survey 
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of the properties of 260 morphologically-selected ETGs within 40 
Mpc (for full details see ICappellari et alj uOllal hereafter Paper 
I). This survey combines integral-field spectroscopy (IPS), opti- 
cal imaging, molecular and atomic gas observa tions and simula- 
tions. As part of this survey lYoung et alj ( 120 lib (hereafter Paper 
IV) have presented an unbiased census of the molecular gas con- 
tent of this sample of nearby early-types, and report that 22% of 
optically-selected, morphologically-classified ETGs have substan- 
tial molecular gas reservoirs (lO^-lO'"' M0 of H2). The brightest 
detections have been followed up with CO-interfeometry (Alatalo 
et al., in prep; see Section [Z2t . Sena et al., in prep report a HI de- 
tection rate of ~40% (for field galaxies with declination above 10° 
drawn from the same sample). These studies are the best available 
complete, unbiased surveys of the ISM in local ETGs, showing that 
a reasonable fraction of these systems have sizable gas reservoirs. 

If galaxies form in a hierarchical manner, then these observa- 
tions pose a challenge to the standard view that early-type galaxies 
join and remain on the red-sequence due to a lack of cold gas. One 
must either demonstrate that it is possible to create a tight red se- 
quence without removing all of the cold ISM, or that galaxies can 
regenerate or acquire cold gas after joining the red sequence. 

The gas we detect in these galaxies could be a remnant, left 
over from the pro genitors that formed the ETG. Stellar popula- 
tion analyses (e.g. iKuntschner et alJuOOa . iKuntschner et al.ll20I(l 
McDermid et al., in prep), however, show that the underlying stel- 
lar populations in these galaxies are old (^ 1 Gyr). When com- 
bin ed with the fact that m ost ETGs show relaxed stellar kinemat- 
ics dEmsellem et alj|201 ll . hearafter Paper III) and evidence from 
deep imaging (Due et al., submitted; Paper IX) this suggests that 
the last major merger that could have formed such systems was at 
least several gigayears ago. Unless the remnant gas can be made 
stable against star-formation this thus suggests that we are seeing 
regenerated or newly accreted gas. 

Galaxies can regenerate a cold ISM through both internal and 
external processes. Stellar evolution models predict that an aver- 
age stellar population returns on the order of half its stellar mass 
to the ISM over a Hubble time 1 Jun^wiert^ Co mbes & Paloua200ll ; 
iLia. Portinari & Carraroll2002l : 



Pozzetti et alj|2007l ; ?). This inter- 



nal gas return is dominated not by cataclysmic events such as su- 
pemovae, but by the long-term stellar mass loss from red giant 
branch, (post-)asymptotic giant branch stars and plan etary nebulae 
l lParriott & Bregmar]|2008l ; iBregman & Parriottll2009l) . The major- 
ity of the ejected mass from giant stars will shock heat and j oin 
the hot gas reservoir of the galaxy jParriott & Bregmanll2008r) . It 
may then be possible for this material to cool, recombine, fall to- 
ward the galaxy centre and eventually become molecular, regen- 
erating the cold gas reservoir. The fraction of gas which is able 
to cool from the hot halo in this way is however unknown. Mass 
loss from planetary nebulae may avoid being shock heated, and 
cool fro m the ionised phase directl y without joining the galaxies 
hot halo dBregman & Parriottll2009h . Simulations have shown that 
this mass loss looses some angular momentum in this process, but 
gene rally it preserves the gener al sense of rotation of the parent star 
(e.g. lMartig & Boumau3l2010h . 

Stellar mass loss must be occurring in all ETGs at all times, 
at a rate depending on the number of stars present in each galaxy 
and its star formation history. One would thus expect many ETGs 
to have detectable molecular gas (assuming that some fraction of 
the hot gas reservoir does cool). Paper IV however found that only 
22% of ETGs have detectable molecular gas, and there is no de- 
pendence of the detection rate and molecular gas mass fraction on 
host luminosity. It is of course possible that gas reservoirs exist in 



the other systems with a more quiescent recent star-formation his- 
tory, resulting in low gas masses from stellar mass loss that are be- 
low our detection threshold. The lack of correlation with the host 
galaxy properties, however, suggests that other mechanisms must 
be destroying or preventing the molecular gas from forming in 78% 
of ETGs, and that stellar mass loss may not be the dominant method 
of acquiring cold gas. 

Galaxies can also regain a cold ISM via external processes, 
such as (major and minor) mergers and/or cold mode accretion from 
the intergalactic medium (IGM). In a minor merger or cold accre- 
tion some fraction of the gas can avoid being shocked and will fall 
to the ce ntre of the galaxy and cool, often producing a c entral disc 
or rings JMazzuca et alj|2006l ; lEUche-Moral et alj[201ol) . In a ma- 
jor merger ejected gas can be re-accreted over the course of cosmic 
time. HI in ETGs has long been thought to come almost excl usively 
from external sources (e.g lKnapp, Turner & Cunniffd 19850 . More 
recently deep HI observations (e.g. lMorganti et al.ll2006h have con- 
firmed that almost all Hl-detected field ETGs show signs of accre- 
tion, indicating an external origin for much of the atomic gas, wit h 
an average HI accretion rate of 0. 1 M© yr~ ^ dOosterloo et alj2010h . 
If the same is true for the molecular and ionised gas has yet to be 
determined. 

Recent work has suggested one way to completely sideste p 
this problem, via 'morphological quenching' dMartig et al.ll2009h . 
In this scenario, red sequence galaxies are created not by removing 
all the cold ISM, but by making some fraction of it stable against 
gravitational collapse. The growth of a stellar spheroid via merg- 
ing or secular processes increases the stellar density, and hence 
the epicyclic frequency and velocity dispersion in t he centre of the 
galaxy. This can cause the Toomre stability criterion dToomrel 19641) 
to exceed unity, implying that the gas is stable against gravitational 
collapse. As the star formation rate would be very low, this stable 
cold ISM would persist for many billions of years, and hence would 
be detectable today in red galaxies. In some ATLAS'^'^ galaxies we 
do indeed see large amounts of HI in a regular disk, bu t no corre- 
sponding young stellar population dMorganti et alj2006h . This pro- 
cess cannot be how every ETG transits to the red sequence, how- 
ever, as such long-lived discs are not observed in all ETGs. 

In summary, determining the dominant source of the cold ISM 
in ETGs is vital in order to understand their formation and evo- 
lution. If stellar mass loss can build up molecular reservoirs, then 
galaxies can transform themselves from spheroidal to disky sys- 
tems over time and in isolation, and perhaps even (if morphological 
quenching does not intervene) evolve back into the blue cloud over 
cosmic timescales (e.g. [Kannappan. Guie & Baker 2009). If how- 
ever mergers are the dominant source of the gas, then star-formation 
episodes are likely to be short-lived and possibly violent, and the 
gas will not accumulate and regenerate. 

Observationally the origin of the gas may be addressed by 
comparing the angular momentum of the ISM with that of the un- 
derlying stellar population. Because of angular momentum con- 
servation, stellar mass loss must produce gas that is kinematically 
aligned with the bulk of the stars which produced it. To first order, 
material from external sources can enter a galaxy with any angular 
momentum, so many such mismatches should be observable today, 
the exact fraction depending on the dynamical timescales for mis- 
aligned material to relax into the equatorial plane. Even in very old 
systems, one might thus expect to see an equal number of galaxies 
with counter and co-rotating gas if external accretion and mergers 
are important (and if there is no preferred accretion direction). 

A number of previous studies have used the misalign- 
ment of the angular momenta of gas and stars to constrain 
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the origin of the ionized gas (e.g. iKannappan & Fabricanj 



S^aetalJ_ 20061) . and the molecular ^as (e.^. [ Young 



2001 



2002 



Young. Bureau & CapDeiiarill2008l : ICrocker et al.ll201lh in early- 
type galaxies. The goal of this paper is to use the ATLAS^'^ com- 
plete volume limited sample of 260 ETGs to constrain the impor- 
tance of externally acquired gas in this way. 

In Section 2 of this paper we describe the observations our 
work is based upon, and present our method for extracting kine- 
matic position angles. In Section 3 we describe the misalignment 
distributions obtained, and we discuss them further in Section 4. 
We conclude in Section 5 and discuss prospects for the future. 

2 Observations and data analysis 

To compare the orientation of the stellar and gaseous angular 
momenta, it is preferable to have full kinematic maps at compa- 
rable angular resolutions. iKrainovic et al.l OOOSh have shown that 
many ETGs have kinematically distinct stellar substructures, that 
may cause confusion in long-slit spectroscopic data. Hence IPS data 
are desirable, where one retrieves full two-dimensional maps of the 
kinematics of the stars and ionised gas. These also allow classifica- 
tion of polar and misaligned structures, instead of confining align- 
ment measures to purely co-rotating and counter-rotating bins (as 
one is forced to do with only long-slit spectroscopy). These 2D 
maps may be combined with mm-wave interferometry, to map the 
cold molecular component of the ISM at similar spatial resolutions. 

The galaxies we use in this work are drawn from the complete, 
volume-limited, morphologically selected ATLAS ^"^ sample (Pa- 
per I) of 260 ETGs brighter than L > 8.2 x 10^ Lq^k, within 
D < 42 Mpc. Paper I summarizes the main data-sets available 
for the survey galaxies. In this paper we combine the ATLAS^'^ 
SAURON optical IPS stellar and ionised gas kinematic maps, CO 
kinematic maps, and HI maps from the Westerbork Synthesis Radio 
Telescope (WSRT). The CO data are a combination of new observa- 
tions from the Combined Array for Research in mm-wave Astron- 
omy (CARMA) and literature data from the Plateau de Bure Inter- 
ferometer (PdBI), Owens Valley Radio Observatory Interferometer 
(OVRO) and the Berkeley-Illinois-Maryland Array (BIMA). These 
observations will be presented in full in Alatalo et al., in prepara- 
tion, but we summarize them briefly below. 

2.1 SAURON IFS data 

SAURON is an IPS built at Lyon Observatory and mounted at the 
Cassegrain focus of the William Herschel Telescope. It is based on 
the TIGER concept (Bacon et al. 1995), using a microlens array to 
sa mple the field of v iew. Details of the instrument can be found 
in iBacon et al.l 1 1200 ih . All galaxies were observed with the low- 
resolution mode of SAURON, covering a field of view of about 
33" X 41" with 0"94 x 0"94 lenslets. Mosaicking was used to 
reach up to one effective radius (the radius encompassing half the 
light) when necessary. 

The SAURON observations for the ATLAS^° galaxies, and 
the extraction of the stellar kinematics are described in detail in Pa- 
per I. In brief, for each tar get, individual datacube s were merged and 
analysed as described in lEmsellem et al.l 1 120041) . ensuring a mini- 
mum signal-to-noise ratio of 40 p er spatial and spectra l pixel using 
the binning scheme developed by ICappellari & Copinl 1 120030 . The 
SAURON stel lar kinematics were derived u sing a penalized pixel 
fitting routine JCappellari & Emsellenj2004r) . providing parametric 
estimates of the line-of-sight velocity distribution for each spaxel. 
During the extract ion of the stellar kinematics, the GANDALF code 
JSarzi et al.ll2006h was used to simultaneously extract the ionised 



gas kinematics. The me an stellar velocity map s for the ATLAS^'^ 
galaxies are presented in lKrainovic et al.l ( 1201 11) . hereafter Paper II. 

2.2 CARMA data 

As part of the ATLAS''^ survey, all galaxies that were detected in 
CO(l-O) (from single pointings of the IRAM-30m telescope at the 
galaxy centre; see Paper IV) with an integrated flux greater than 19 
Jy km s~^ that do not have interferometric data already availa ble i n 
the literature have been observed with CARMA (iBock et aUi2006h . 
Full details of this interferometric survey can be found in Alatalo et 
al. (in preparation), but we summarize the observations here. 

Observations of the detected sample galaxies were taken be- 
tween 2008 and 2010, mainly in the D-array configuration, pro- 
viding a spatial resolution of 4-5". CO(l-O) was observed using 
narrow-band correlator configurations, providing at least 3 raw 
channels per 10 km s~^ binned channel whilst ensuring adequate 
velocity coverage for all galaxies. 

Bright quasars were used to calibrate the antenna-based gains 
and for passband calibration. The data were calibrated and im- 
aged using the Multichannel Image Recons truction, Image Analysis 
and D isplay (MIRIAD) software package dSault. Teuben & Wrighll 
1199 51) . Total fluxes and velocity widths were compared with the 
IRAM 30m single-dish observations to ensure that large propor- 
tions of the fluxes were not being resolved out. 

A total of 27 galaxies included in this work have been 
observed and all were detected at greater than 3ct with CARMA 
(IC676, IC719, IC1024, NGC 1222, NGC2697, NGC2764, 
NGC2824, NGC 3182, NGC 3607, NGC 3619, NGC 3626, 
NGC 3665, NGC 41 19, NGC 4324, NGC 4429, NGC 4435, 
NGC 4694, NGC 47 10, NGC 4753, NGC 5379, NGC 5866, 
NGC 6014, NGC 7465, PGC 029321, PGC058114, UGC 06176 
and UGC 09519). We also include here galaxies for which 
CO(l-O) maps are alre ady available from th e literature, mostly 
from SAURON survey dde Zeeuw et al.l l2002h follow-ups. These 
are NGC 0524, NGC 2685, NGC 2768, NGC 3032, NGC 3489, 
NGC 4 150, NGC 4459, NGC 4476, NGC 4477, NGC 4526 and 
NGC 4550. Full references for the data used in this work are listed 
in Table [T] This makes for a total of 38 sample galaxies that have 
CO(l-O) interferometric data. 

This sample of mapped galaxies are the brightest ~2/3 of the 
CO detections from Paper IV. The mapped galaxies are thus bi- 
ased towards higher molecular gas masses, but are statistically in- 
distinguishable (with a Mann-Whitney U (MW-U) or Kolmogorov- 
Smirnov test) from the full sample of CO detected systems in 
terms of their molecular mass fractions, local environment, and host 
galaxy properties. Th e mapped galaxie s also fall on the same CO 
Tully-Fisher relation l IDavis et alJl201lL hearafter Paper V). Other 
than molecular gas mass, the only variable in which our sample 
may have some bias is in cluster membership. Due to the fact that 
the Virgo cluster is relatively nearby (and hence the flux from a 
given amount of molecular gas is higher) we have mapped a greater 
percentage of cluster members than field galaxies. We do not how- 
ever expect this to affect our conclusions. Overall we consider that 
the CO sample discussed from this point should be reasonably free 
of biases, however it is worth remembering that we do not full sam- 
ple the parameter space at smaller molecular gas masses. 

In this work we make use of the CO velocity fields which 
will be presented in Alatalo et al, in prep. These mean veloc- 
ity maps were produced by the masking method: each fully cal- 
ibrated and cleaned image cube was smoothed along both the 
spatial and velocity axes, and the smoothed cube was clipped at 
~2.5 times the rms noise in a channel. The clipped version of 
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the smoothed cube was then used as a mask to define a three- 
dimensional volume in the original, unsmoothed cleaned cube, over 
which we calculate the first velocity moment (see lRegan et al.l200lK 
[Young. Bureau & CaoDellari 2008). An identical method was used 
inB'ouns. Bureau & Ca ppellaril ( l2008l) : ICrocker et al.l l l2008ll2009h 
and .Crocker et al.. C20I1|) . and we refer readers to these works for 
examples of the data quality achieved. 

2.3 Measuring kinematic position angles 

The apparent (projected) direction of the angular momentum vec- 
tor of any tracer can be found by estimating its kinematic position 
angle 0, which lies normal to the angular momentum vector. The 
kinematic position angle (PA) is defined as the counter-clockwise 
angle between north and a line which bisects the tracer's velocity 
field, measured on the receding side. 

One can then define the kinematic misalignment angle ip as the 
difference between the kinematic PA of two galactic components. In 
this paper, we will consider misalignments between the molecular 
gas, warm ionised gas and stellar angular momentum; 



V^mol — star ^ | fpmolccular gas fpstar | , 
'Pmol — ion ^ | v^molccular gas v'ionised gas | , 
V^ion — star =^ | ff^ioniscd gas ff^star | , 



(1) 

(2) 

(3) 



These quantities are defined to lie in the range 0-180°. 

In this work, the kinematic PA of the stars is taken from Pa- 
per II (Table 1, column 4 in Appendix C, and available online at 
|http ://www.purl.org/atlas3d I. The molecular and ionised gas kine- 
matic PAs are presented here for the first time, and are measured 
using the same m ethod (the FIT_K INEMATIC JA routinqa described 
in Appendix C of lKrainovic et alj|2006) . 

In this paper we were able to measure molecular gas kinematic 
position angles for all the galaxies with mapped molecular gas. The 
values of the kinematic PA for these galaxies are listed in Table [T] 
We give two examples of CO velocity fields overlaid with their de- 
rived kinematic PA and the associated error in Figure[T] Ionised gas 
kinematic position angles were measured only on galaxies where 
the emission appeared coherent, with a regular velocity field, and 
was not very confined to the very central regions (as this is in- 
dicative of templ ate mismatch within the GANDALF package; see 
ISarzi et alJuOOq) . These criteria were satisfied in 132 of the 260 
sample galaxies (51%). Where the ionised gas is only detected over 
a small region, we have re-run the FIT_KINEMATIC_PA code using 
only this region to reduce the effect of low signal to noise ratio bins 
on our PAs determinations. Table [T] contains the ionised gas kine- 
matic PA for all CO mapped galaxies in this work. Tables [3] and 
and |4] contain the kinematic misalignment between the ionised gas 
and stellar components for all fast-rotating galaxies and slow rotat- 
ing galaxies respectively, which were not mapped in CO, but which 
have a significant enough ionised gas detection to allow a kinematic 
PA measurement. 

HI data was available for 29 fast-rotating galaxies that over- 
lap with the CO and ionised gas sample considered in this paper. 
The large difference in angular scale probed make comparisons be- 
tween the extended HI and centrally concentrated molecular and 
ionized gas challenging, but potentially powerful. We do not extract 
kinematic position angles for the HI, but discuss the general trends 
present from visual inspection of the HI velocity fields in Section 
13.41 More detailed analysis of the HI velocity fields and intensity 
maps will take place in a future paper in this series. 



§ Available from http://purl.org/cappellari/idl 
© 2009 RAS, MNRAS OOO.mfTTl 



In real galaxies, there can be kinematic substructures (e.g. 
bars, kinematically decoupled cores) that may or may not share the 
same kinematic PA. In these cases, the fit Jcinematicpa routine re- 
turns the value applicable to the bulk of the material (as it is area, 
rather than luminosity weighted). Such substructures can, however, 
skew the measured kinematic PA, in some cases by tens of degrees. 
When a bar is present. Paper II has shown that the global kinematic 
PA of the stars is uncertain with a standard deviation of ~10°. In 
general however the stellar kinematic position angle is a good es- 
timate of the PA of the line of nodes even when a bar exists. Bars 
would be expected to have a stronger influence on the gaseous com- 
ponents than on the stars. 

It is worth remembering that we always compare the kinematic 
PAs of the molecular gas to that of the main body of the stars. Some 
stellar velocity fields show a kinematically decoupled core (KDC), 
a stellar subcomponent with misaligned rotation from the rest of the 
stellar body. These KDCs generally extend only over a small central 
region however, and hence do not significantly affect our results. 

The mm interferometric data have a poorer angular resolu- 
tion than the optical data (~4-5" vs 1-2"). In order to test how 
this affects the determination of the kinematic PA, we performed a 
Monte-Carlo simulation. A synthetic velocity field with a kn own PA 
was cr eated using an empirical galaxy rotation curve from lRoscod 
1 119991) . We then added Gaussian noise to simulate the average sig- 
nal to noise ratio of our data. We varied the inclination, and the 
angular extent of the emitting material while keeping a fixed beam 
size (which is identical to changing the angular resolution) using a 
fixed spatial sampling of 3 pixels per beam, and measured the kine- 
matic PA using the fit Jcineinaticpa routine. Changing the inclina- 
tion from edge-on to face-on had little effect of the determination 
on the PA, apart from very face on inclinations where the veloc- 
ity range becomes comparable to the channel width. This happens 
only at inclinations less than 2.8 degrees for a typical galaxy with a 
maximum velocity of 200 km s~^, and 10 km s~^ channel width. 
Varying the source size from 50 to 2.5 beams across (or equiva- 
lently 1" to 20" angular resolution) caused a scatter of ~10° in the 
determined kinematic PA, with no systematic offset. The true PA 
was always within the 3a error bars returned by the routine. 

We also used this simulation to test the effect of signal to noise 
on our PA determinations. The average RMS deviation from the 
bi-(anti)symmetric velocity map created by the fitjcinematicpa is 
«15 km s^^ for the real CARMA data used in this work (with a 
maximum of 35 km s~^). Varying the simulated RMS noise from 5 
km s~^ to 50 km s~^ scattered the derived kinematic PAs by ~10°. 
Once again this induced no systematic offset, and the true PA was 
always within the 3a error bars returned by the routine. 

The largest error on the observed molecular gas PA (as esti- 
mated by the fit Jcinematic. pa routine) is 30°, and the average (Icr) 
error on the molecular gas PAs is 7°, consistent with the predicted 
error from our simulation. In this paper we choose to consider kine- 
matic misalignments between the gas and stars of > 30° as being 
significantly misaligned. This conservative estimate ensures that we 
can set a robust lower limit on the importance of externally accreted 
material. 

3 Results 

3.1 Molecular gas and stars 

A direct comparison between the kinematic PAs of the molecular 
gas and the main body of the stars reveals two populations (Figure 
|2ll. We find 24±7% (9 out of 38 galaxies) have their molecular gas 
misaligned ('l/ > 30°) with respect to the stars, setting a lower 
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Figure 1. Examples of the data for one of the galaxies with the best constrained molecular gas PA (NGC 0524, top), and the galaxy with the w orst constrained 
molec ular gas PA (NGC2768, bottom). Left: CO velocity field (colours). Mid dle: Stellar velocit y fields from the SAURON IFU (data from lEmsellem et alj 
|2004J). Right: H/3 emission line kinematics from the SAURON IFU (data from lSarzi et al.ll2006l) . These are overlaid with the CO integrated intensity (black 
contours) and the best fitting molecular gas/st ellar kinematic PA, and its associate d error as calculated in this paper (solid black line, and dashed black lines 
respectively). The CO moments are taken from lCrocker et al.l2008l (NGC2768) and lCrocker et al.l201 ll (NGC0524). The SAURON IFU coverage of NGC2768 
extends beyond the range of the plots shown above. 



limit on the number of galaxies where the gas is likely to have an 
external origin. The remaining 76±7% (29 out of 38 galaxies) have 
their molecular gas aligned with the stars (^ < 30°) with the stars, 
consistent with an internal origin. The uncertainties quoted above 
include only the formal binomial counting errors. 

Several individual galaxies have a KDC coincident with the 
molecular gas. NGC 4476 is a pathological example, as it has a 
large stellar KDC, kinematically aligned with (and likely formed 
from) the molecular gas, while the stars have little sense of ro- 
tation outside this region within the SAURON IFS field of view. 
PPAK IFS stellar kinematics suggest that the stellar velocity field 
may become inverted outside the SAURON field of view (Crocker 
et al., in prep), hence the molecular gas could be considered to be 
counter-rotating with respect to the stellar body at large radii. For 
consistency, however, we use only the SAURON stellar kinematics 
in this work (probing roughly one effective radius), and hence la- 
bel NGC 4476's molecular gas as co-rotating. It is worth bearing in 
mind, however, that the stellar kinematics in the outer parts of all 
of these galaxies could be different from the kinematics in the cen- 
tral regions. Discounting the few galaxies such as this, where the 
origin of the gas is hard to judge would not significantly alter our 
conclusions. 

3.2 Molecular and ionised gas 

The analysis presented in Section [3711 can be repeated, comparing 
the apparent angular momenta of the molecular and ionised gas. 
This results in Figure [3] The molecular gas has a kinematic mis- 
alignment angle with the ionised component of less than 40° in 
every galaxy (and generally much better than this), and within the 
errors quoted all galaxies are consistent with having a kinematic 
misalignment angle of less than 10.5°. The overall distribution is 
consistent with a Gaussian error distribution with centre at zero and 
a standard deviation of 15°, as shown in Figure [3] 
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Figure 2. Histogram showing the kinematic misalignment angle between 
the molecular gas and the stars for all of the galaxies listed in Tabled The 
dashed line overplotted is a normalized Gaussian distribution with its center 
at zero and a standard deviation of 7°, showing the expected scatter 



The galaxies which have a misalignement between 15° 
and 40° are IC0676, NGC 1222, NGC 3032, NGC 4150 and 
PGC 029321. These small misalignments can always be under- 
stood as arising from small-scale kinematic substructures within 
the galaxies affecting the ionised component, that are usually un- 
resolved at the resolution of our CO interferometry. IC 0676 has its 
ionised gas kinematic PA affected by a bar. NGC 1222 is a highly 
disturbed star-bursting sys tem, that may be in the process of un der- 
going a three-way merger l IBeck. Turner & Kloostermanll2007h . As 
discussed above, NGC 3032 an d NGC 4150 have KDCs that may 
be affecting the gas properties l IMcDermid et al.ll2006l ; ISarzi et alj 
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Table 1. Kinematic parameters of tlie ATLAS^^ early-type CO-rich galaxies used in this paper 



Name 


(/•mol 


A0mol 


0ion 


A</.ion 


V'mol-star 


AV'mol-star 


V^mol — ion 


AV^mol-ion 


V'ion — star 


AV^ion-atar 


Ref. 




(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


IC()676 


16.5 


10.0 


348.5 


26.9 


2.0 


15.8 


28.5 


29.6 


30.0 


29.6 




IC0719 


229.0 


3.5 


232.5 


2.5 


177.0 


26.2 


3.5 


26.1 


173.5 


26.1 




IC1024 


24.5 


12.8 


31.5 


8.2 


5.0 


16.7 


7.0 


13.6 


2.0 


13.6 




NGC0524 


36.5 


1.8 


41.5 


21.9 


4.0 


2.7 


5.0 


22.0 


1.0 


22.0 




NGC1222 


33.0 


1.5 


54.5 


21.5 


10.0 


9.4 


21.5 


23.4 


11.5 


23.4 




NGC2685 


105.0 


10.0 


109.5 


13.8 


68.5 


10.3 


4.5 


14.0 


73.0 


14.0 




NGC2697 


301.5 


1.8 


300.0 


3.1 


0.0 


3.6 


1.5 


4.4 


1.5 


4.4 




NGC2764 


202.5 


1.8 


189.5 


5.2 


6.5 


7.0 


13.0 


8.6 


6.5 


8.6 




NGC2768 


180.0 


30.0 


187.5 


6.1 


87.5 


30.2 


7.5 


7.0 


95.0 


7.0 


4 


NGC2824 


161.5 


1.8 


159.5 


7.5 


2.0 


3.3 


2.0 


8.0 


0.0 


8.0 




NGC3032 


92.5 


10.0 


62.5 


19.1 


179.0 


14.9 


30.0 


22.0 


151.0 


22.0 




NGC3182 


331.5 


7.2 


316.5 


6.0 


11.0 


10.0 


15.0 


9.2 


4.0 


9.2 




NGC3489 


67.5 


1.8 


78.5 


5.3 


5.0 


3.3 


11.0 


6.0 


6.0 


6.0 




NGC3607 


302.5 


2.5 


302.5 


0.5 


1.0 


2.9 


0.0 


1.6 


1.0 


1.6 




NGC3619 


74.5 


4.2 


76.0 


0.5 


22.5 


5.2 


1.5 


3.0 


24.0 


3.0 




NGC3626 


169.5 


1.8 


165.0 


3.2 


170.0 


1.9 


4.5 


3.3 


174.5 


3.3 




NGC3665 


219.5 


2.0 


209.5 


2.2 


14.0 


2.8 


10.0 


3.0 


4.0 


3.0 




NGC4119 


296.0 


12.0 


292.5 


18.6 


4.5 


13.6 


3.5 


19.7 


1.0 


19.7 




NGC4150 


146.0 


10.0 


168.5 


10.1 


1.5 


11.9 


22.5 


12.0 


21.0 


12.0 




NGC4324 


232.0 


1.8 


239.0 


6.8 


6.0 


5.5 


7.0 


8.6 


1.0 


8.6 




NGC4429 


82.0 


2.0 


92.5 


2.0 


4.5 


3.2 


10.5 


3.2 


6.0 


3.2 




NGC4435 


201.0 


1.8 


198.6 


3.1 


8.5 


2.5 


2.4 


3.6 


6.1 


3.6 




NGC4459 


269.0 


5.0 


281.5 


1.5 


11.5 


5.6 


12.5 


2.9 


1.0 


2.9 


3 


NGC4476 


208.0 


2.5 


219.0 


7.5 


1.5 


11.8 


11.0 


13.7 


12.5 


13.7 


1 


NGC4477 


227.0 


3.0 


226.5 


1.6 


25.5 


6.0 


0.5 


5.5 


26.0 


5.5 


6 


NGC4526 


288.0 


3.0 


294.5 


4.7 


0.5 


3.5 


6.5 


5.0 


6.0 


5.0 


3 


NGC4550 


355.0 


10.0 


358.5 


9.7 


3.5 


10.3 


3.5 


10.0 


0.0 


10.0 


5 


NGC4694 


155.5 


27.5 


167.0 


45.0 


169.0 


33.6 


11.5 


48.9 


157.5 


48.9 




NGC4710 


207.0 


10.0 


208.5 


2.2 


0.5 


10.7 


1.5 


4.4 


1.0 


4.4 




NGC4753 


93.0 


2.8 


85.5 


1.5 


4.5 


3.8 


7.5 


2.9 


3.0 


2.9 




NGC5379 


66.0 


20.0 


77.0 


5.0 


5.0 


22.4 


11.0 


11.2 


16.0 


11.2 




NGC5866 


127.0 


2.5 


123.0 


4.5 


2.0 


3.9 


4.0 


5.1 


2.0 


5.4 




NGC6014 


139.5 


6.2 


151.5 


10.2 


7.5 


10.5 


12.0 


13.3 


4.5 


13.3 




NGC7465 


106.0 


3.8 


109.5 


9.5 


60.5 


29.2 


3.5 


30.5 


57.0 


30.5 




PGC029321 


76.0 


27.5 


105.5 


27.2 


19.5 


46.1 


29.5 


46.0 


49.0 


46.0 




PGC058114 


94.5 


7.0 


93.1 


3.1 


152.5 


12.2 


1.4 


10.5 


153.9 


10.5 




UGC06176 


201.0 


8.5 


214.0 


8.5 


0.5 


30.2 


13.0 


30.2 


13.5 


30.2 




UGC09519 


177.5 


5.5 


172.0 


4.5 


72.0 


8.0 


5.5 


7.3 


77.5 


7.3 





Notes: Columns 2-5 list the kinematic PAs for the molecular gas and ionised gas (even columns), and their respective errors (odd columns). Columns 6-11 list 
the misalignments between the molecular gas, ionised gas and stars (even columns) and and their respective errors (odd columns). The errors quoted on the 
kinematic PAs are those calculated by the FIT_KINEMATIC_PA routine. The kinematic PA of the main body of the stars for each galaxy is taken directly from 
Paper II (Table 1 in Appendi x C, Column 4; available online at http:// ww w.purl.o rg/atlas3d I. Colu mn 12 lists the data re fere nce: (1) Alatalo et al. , in prep (2) 
ISchinnerer & Scovillel ^002l) . (3) lYoung. Bureau & Cappellaril J2008l) ,T4) ICrocker et al, l^QOj), (5) ICrocker eTZI )2009l) . (6) ICrocker eTaJ] )201 ih . 



l200a) . The ionised and molecular material in PGC 029321 has a 
small angular extent, and as we showed above this can increase the 
error in the kinematic PA determination. 

The uncertainty on ionised gas PA measurements is larger than 
the average 10° we have estimated for the molecular gas PAs. As 
discussed above, this is because the ionised gas is often patchy, has 
a smaller angular extent than the molecular gas, and can be less re- 
laxed. Despite this increased uncertainty, and greater intrinsic vari- 
ance due to small substructures in the galaxies, our conservative 
threshold for misalignment (\I' > 30°) ensures that we still set a 
strong lower limit on the importance of externally accreted material. 
Despite the increased uncertainty, we only show plots of the ionised 
gas misalignment from here on, taking advantage of the increased 
number of galaxies available to improve our statistics. When one 
examines the molecular gas kinematic misalignments all results re- 
main unchanged. 

One galaxy (NGC 3032) has a molecular to ionized gas kine- 



matic misalignment > 30° cut-off we have taken to denote the 
boundary between aligned and misaligned gas. As the distribu- 
tion of misalignments closely follows the gaussian error distribu- 
tion however, we suggest that this is due to observational scatter. 
We therefore conclude that the ionised and molecular material are 
always kinematically aligned with each other, even when they are 
kinematically misaligned with respect to the stars, strongly suggest- 
ing they share a common o rigin (as previously discussed in a sub- 
sample of these galaxies bv lCrocker et al.ll201lh . This is discussed 
further in Section|4T| 

Assuming 30° rather than 40° for the cut-off between aligned 
and misaligned gas can at most bias our results by ~4% (this is the 
area under the error distribution at misalignments greater than 30°, 
normalized by the total area under this error distribution curve.) 
For the typical number of galaxies included in our analysis this can 
result in the misclassification of 1-2 galaxies. 
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Figure 3. Histogram sliowing tlie kinematic misalignment angle between 
the molecular gas and the ionised gas for all the galaxies listed in Tabled 
The dashed line overplotted is a normalized Gaussian distribution with its 
center at zero and a standard deviation of 15°, showing the expected scatter. 



Figure 4. Histogram showing the kinematic misalignment angle between 
the ionised gas and the stars for all the slow-rotating galaxies with measur- 
able kinematic misalignments. These systems are listed in Table|4] 



3.3 Ionised gas and stars 

As the molecular and ionised gas are always aligned and are likely 
to share a common origin, investigating the kinematic misalignment 
between the ionised gas and the stars can therefore also shed light 
on the origin of both of these gas phases in ETGs. 

The ATLAS'^^ ionised gas detection rate (above an integrated 
equivalent width of 0.02 A in either H/3 or [OIII]) is Ri73±3% in the 
field, and «47±6% in the Virgo cluster. These figures are compa- 
rable with the previous SAURON survey ionised gas de tection rate 
of 86 ±6% in the field and 55±12% in the Virgo cluster JSarzi et alj 
120061) . We are hence able to repeat the analysis in Section lSTl with 
the benefit of increased number statistics. 

Table [3] contains the kinematic misalignment between the 
ionised gas and stars for all fast-rotating galaxies with a significant 
enough ionised gas detection to allow a kinematic PA measurement. 

3.3.1 Slow-rotators 

Slow-rotators as a class are round, massive, dispersion dominated, 
mildly triaxial ETGs that generally have old stellar populations, 
and large KDCs ( Paper III). Ionised gas is detected in slow-rotators 
l lSarzi et alj|2006f) . while molecular gas almost never is (Paper IV). 
Ionised gas is hence the best tracer to help us understand the origin 
of the gas in the central parts of these systems. The stellar kinematic 
PA can be poorly defined in a slow-rotator, which have little or no 
coherent rotation by definition, and this translates to large uncer- 
tainties in the misalignments. Despite this it is possible to estimate 
the kinematic misalignment for 22 of the 36 slow rotating galaxies 
in the ATLAS'^'^ sample. The derived misalignments are listed in 
TableU 

When slow rotators are considered by themselves (Figure 0, 
they display a flat distribution of misalignments between the ionised 
gas and the stars, suggesting the dominant source of (io nised) gas is 
extern al (as was seen previously with fewer galaxies bv lSarzi et alj 
120071) . A MW-U test finds no statistically significant evidence that 
the observed kinematic misalignments (in Figure |4ll are not drawn 
from a uniform underlying parent distribution. 

3.3.2 Fast-rotators 

Figure |5] shows that 36±5% of fast-rotating galaxies (40 out of 
111) have their ionised gas kinematically misaligned from the stars 
($ion-star > 30°), consisteut with the percentage of misaligned 
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Figure 5. Histogram showing the kinematic misaUgnment angle between 
the ionised gas and the stars for all the fast-rotating galaxies listed in Table 
[3] The dashed line overplotted is a normalized Gaussian distribution with its 
center at zero and a standard deviation of 15°, showing the expected scatter 



galaxies found in molecular gas (but with a smaller uncertainty). 
Where both the ionised gas and molecular gas are detected, the mea- 
sured kinematic misalignments from the stars agree well with each 
other, with a la scatter around the one-to-one relation of 10°. A 
MW-U test finds no statistically significant evidence that the molec- 
ular gas and ionised gas misalignments from the stars (Figures[2]and 
|5} are not drawn from the same underlying parent distribution. 

3.4 Comparison to HI 

As part of the ATLAS '^'^ survey we collected HI data for all 166 
galaxies above declination -1-10° and further than 15 arcmin from 
Virgo A. Most of these galaxies were observed with the Westerbork 
Synthesi s Radio Telescope (WSRT). Observations are d escribed in 
detail in iMorganti et al.l ||2006|): lOosterloo et alj | |2010|) and Serra 
et al., in prep (see ISerra et alJl2009L for a summary). For about 
20 galaxies inside the Virg o cluster we use Arecibo data taken as 
part of the Alfalfa survey jdi Serego Alighieri et alj|2007h instead 
of WSRT data. We use Alfalfa data to derive an upper limit on the 
HI mass which is comparable to what we obtain with our WSRT 
observations. 
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This ATLAS HI sample has 29 HI detections of fast-rotators 
that overlap with the molecular and ionised gas sample with mea- 
surable misalignments we consider in this work (13 that overlap 
with the interferometric CO sample, and a further 16 that overlap 
with the ionised gas sample). Despite the large difference in angular 
resolution between the HI data and the other tracers, in most cases it 
is possible to estimate a rough kinematic PA of the atomic gas in the 
region where it overlaps with the CO and ionised gas. We find that 
the atomic and molecular material always appear to be kinemati- 
cally aligned with one another, suggesting that the y share a com- 
mon o rigin. This is in agreement with the results of lMorganti et alj 
<200q) . who find that the HI and ionised gas are always aligned in 
the SAURON sample (which is a subset of the ATLAS^° sample). 

As we have established that the ionised, molecular and 
atomic gas are all likely to be part of the same structure, it 
is interesting to use the HI data to investigate the connection 
with HI seen at large radii. Of the Hl-detected systems that 
have (molecular or ionised) gas kinematically misaligned from 
the stars, 86% (12/14) have extended HI structures, clearly sug- 
gesting a source for the misaligned ma te rial. T hese systems 
are discussed in detail i n iLi & Seagiiis^ (1 1994) (NGC7465) , 



HI detections in Virgo (e.g. Giovanar di. Rrum m & SalDetenll983l : 



IVerheiien & ZwaanI ( 1200 ih (NGC 4026). iMorganti et"!!. 
(NGC 1023, NGC 2685, NGC 2768, NGC 7332). IDuc et alj l 



200€) 



120071) 



(NGC 4694). lOosterloo et al] ( |2010|) (NGC 3032. NGC 4262), and 
Serra et al., in prep (NGC 0680, NGC 3626, UGC 09519). 

Furthermore, 38% (6/16) of the Hl-detected galaxies that have 
gas kinematically aligned with the stars also have kinematically dis- 
turbed HI at larger radii, with clear suggestio ns of ongoing accretion 
from external sources. These are disc ussed in lMorganti et alj ( l2006h 
(NGC 4150). lOosterloo et al] ilOld) (NGC 3384, NGC 3489) and 
Serra et al., in prep (NGC 2778, NGC 3457, PGC 028887). 

3.5 Environmental effects 

3.5.1 Cluster environment 

The ATLAS'^° sample includes a range of galaxy densities, from 
the Virgo cluster to sparse fi eld environments. Galaxies in clus- 
ters are often HI deficient (e.g. ldi Serego Alighieri et alj2007r) . and 
it is thought that the atomic material is removed by ram-pressure 
stripping and interactions as the galaxies fall into the cluster po- 
tential and plow through t h e hot IC M (e.g. Gunn & Go"t^ 1 19721 : 
iGiovanelli & Havnedl 19831 ). iKennev & Yound ( Il986h have shown 
that spiral galaxies in Virgo are not CO deficient, even when they 
are HI deficient. This is likely to be because molecular gas is denser 
and resides deeper in a galaxy potential well, where it will not be as 
affected by infall into a cluster. 

To first order, the merger rate is enhanced as galaxies fall into 
clusters, but once the galaxies are virialized the chance of a merger 
drops off considerably (the large velocity dispersion within the clus - 
ter makes mergers reasonably rare; e.g. Ivan Dokkum et al.lll999l) . 
Cluster (and group) centre galaxies are the only exceptions, which 
are likely to have many merging events, and are typically slow ro- 
tators (see Paper VII). Paper IV shows that ATLAS^° Virgo fast- 
rotating ETGs are likely to be virialized within the cluster (the line- 
of-sight velocities are centrally peaked and consistent with being 
relaxed), and hence their last merger should have occurred several 
crossing times ago (a few gigayears). 

Stripped or free-floating gas is likely to be destroyed and shock 
heated to the cluster virial temperature as it is captured by the clus- 
ter, joining the cluster's hot X-ray halo. The external channels by 
which cold gas can be acquired are therefore likely to be shut-off 
once a galaxy is embedded in a cluster potential, whereas internal 
stellar mass loss is not impeded. This is consistent with the lack of 



Idi Serego Alighieri et aUl2007l : lOosterloo et al.ll20ld and Serra et 
al., in prep) and with the steep decrease in the frequency of spirals, 
and increase in the frequency of fast-r otating ETGs in the Virgo 
cluster core (see lCappellari et alJuOl Ibl hereafter Paper VII) 

Despite this, as for spirals, the ETG CO detection rate inside 
the Virgo cluster is much higher than the HI detection rate, and is 
almost identical to that found in the field (Paper IV). The molecu- 
lar gas mass fractions (Mhs/Lk) of cluster and field galaxies also 
appear to be similar. Paper IV does however show some evidence 
that the H2-richest ETGs are all low mass, and in poor environ- 
ments. It is hence interesting to investigate these high- and low- 
density galaxy populations separately, to see if the origin of the gas 
is affected by the environment. 

When our sample of galaxies is divided into Virgo (defined 
in Paper I as all the galaxies within a sphere of radius 3.5 Mpc 
centered on the cluster centrq^ll) and field (everything else), a clear 
dichotomy in the kinematic misalignment distribution is apparent. 

In Figure m 42±5% (38 out of 91) of the fast-rotating field 
ETGs have kinematically misaligned ionised gas. However, in 
the Virgo cluster only 10±6% (2 out of 20) of the fast-rotating 
galaxies have ionised gas kinematically misaligned for the stars 
(*ion-star < 30°). If this distribution of alignments in the clus- 
ter was the same as in the field one would expect ~8 to 10 galaxies 
to show misalignments. 

A MW-U test gives a probability that the Virgo and field galax- 
ies are randomly drawn from the same parent galaxy population as 
0.3% (~3o- significance). The same result is found when one con- 
siders the kinematic misalignment between the molecular gas and 
the stars, but with lower number statistics (shown as the hatched 
area in Figure |6j. 

The two kinematically misaligned Virgo galaxies (NGC 4262 
and NGC 4694) both lie on the outskirts of the cluster and are under- 
going mergers/interactions as they fall into the cluster. NGC 4262 
may have recently undergone a close encounter with a nearby 
galaxy (see IVollmer. Huchtmeie r & van Driel' 20051 ). It has a ring 
of HI surrounding the galaxy (Krumm, van Driel & van WoerdenI 
Il985n . This system is strongly barred, and the HI shows elliptical 
orbits which may be the cause of the observed misalignments. This 
ring may be the source of the ionised gas as the kinematic position 
angles appear to be similar. No CO emission was detected from 
NGC 4262 (Paper IV). NGC 4694 is undergoing a gas-rich major 
merger in the cluster outskirts. The CO and HI in this galaxy are 
in non-equilibri um structures of fset from the stellar body, towards 
a dwarf galaxy JDuc et alj2007l : Alatalo et al, in prep.). 

3.5.2 Group environments 

There is observational evidence that preprocessing (through 
mergers and tidal disruption) in group environments may be 
important in the formation of ETGs i Zabludoff&_Mulchae^ 



1998|jKodama et al.l200ll : lHelsdon & Ponmadl2003l : lKautsch et'al 



20081) . Paper VII also shows that local environment density seems 
to be the most important factor in determining the fraction of spiral 
galaxies and fast-rotating ETGs, outside of the Virgo cluster core. 
However, simulati ons show that galax ies can enter clusters directly 
from the field (e.g. lSerrier et al.ll2009l) . avoiding any preprocessing 
in groups, and that therefore any transformation mechanism should 
also be able to occur in the cluster environment. 

It is possible with the data we possess to examine if it is the 



^ Cluster membership is listed in Paper I, Table 5, Column 6, and available 
online at|http://www. purl.org/atlas3d | 
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Figure 6. Top: Histogram sliowing tiie idnematic misalignment angle between the ionised gas and the stars for all fast-rotator galaxies in Virgo (defined as lying 
within a sphere of 3.5 Mpc centered on the cluster centre). The hatched area indicates the number of galaxies in each bin that were also mapped in molecular 
gas. The two outlying galaxies with ^t >30 are NGC4262 and NGC4694. Bottom: As above, but for all fast-rotator galaxies which are not in Virgo, including 
galaxies in the field and small groups. 



local density or the global environment that drive the observed 
kinematic misalignment correlations. To do this, we utilize a lu- 
minosity surface density estimator (I3) based on a redshift cylin- 
der with a depth of 600 km s~^, and an angular size adapted to 
include the 3rd nearest neighbor, as presented in ICappellari et al.l 
( 1201 Ibf) (Paper VII; Table 2, Column 7 and available online at 
|http://ww w.purl.org/atlas3d I. This estimator should pick out locally 
dense environments, such as groups, as well as substructures within 
the Virgo cluster. We plot this density indicator against the kine- 
matic misalignment of the ionised gas and stars in Figure|7] 

In the densest environments, above a critical local luminosity 
surface density of « 10^^'^ Lq Mpc~^, all fast-rotators in both the 
field and cluster have ionised gas kinematically aligned with the 
stars, while a wide range of misalignments exist at lower densities. 
The two misaligned systems in the Virgo cluster are below this lu- 
minosity surface density (they are thought to be recently accreted 
systems, still undergoing mergers in the cluster outskirts; see Sec- 
tion l3.5.ll . This result does not change if one plots number surface 
density (E3; see Paper VII, Table 2, Column 4) rather than lumi- 
nosity surface density, with the critical number density being ~ 20 
Mpc ~ ■^ . A MW-U test gives a (a posteriori) probability of ~ 1 % that 
the kinematically aligned galaxies above (and all the galaxies be- 
low) this critical density are randomly drawn from the same parent 
population. If one considers the larger scale Jio and Eio environ- 
ment indicators from Paper VII the cluster/field result discussed in 
Section [3.5.1l is recovered. 

The field galaxies that exist above this critical density 
are NGC0524, NGC0680, NGC3379, NGC3412, NGC5353, 
NGC 5355, NGC 5379 and PGC 042549, all of which are in known 
galaxy groups. All but one of these systems (which were observed 



in our HI survey) are HI poor (NGC 0680 is the exception, it is 
a major merger remnant, and has kinematically disturbed HI ex- 
tending between several of its group members; see Paper IX and 
Serra et al., in prep). A t least two of thes e groups have extended 
X-ray halos (NGC 524. IJomer et alj|2o"ool ; NGC 5353/NGC 5355, 
lO'Sullivan. Forbes & Ponmadl200lh . 

3.6 Mass dependence 

Whenever one sees an effect which appears correlated with envi- 
ronment, it is important to check if the observed relation is in fact 
driven by galaxy mass. The mass of a galaxy can affect the gas con- 
tent in various ways. A deeper potential can allow a galaxy to retain 
gas that would be lost due to supernovae and galactic winds in a 
lower mass system. Massive ETGs can also have X-ray halos, which 
can prevent external gas (especially in cold flows) from entering 
the galaxy. Figure [8] shows the dependence of the ionised gas-star 
kinematic misalignment of fast-rotators on absolute iCs-band mag- 
nitude, an observable proxy for stellar mass. Distances to the galax- 
ies are taken from Paper I. One can see a clear lack of misaligned 
fast-rotating galaxies at Mk ^ — 24. The only kinematically mis- 
aligned galaxy brighter than this threshold is NGC 2768, which is 
a large elliptical galaxy in a relatively isolated en vironment, that 
appears to be accreting fr om a nearby cloud of HI JMorganti et alj 
l2006l ; ICrocker et al.ll2008h . A MW-U test gives a «2% chance that 
the galaxies with Mk ;$ — 24 are drawn from the same distribution 
of misalignment angles as the fainter galaxies. 

We re-plot Figure |7] with the galaxies with Mx < — 24 mag 
highlighted in green (Figure|9). Interestingly, the fast-rotators above 
the proposed critical luminosity are not all the same systems iden- 
tified previously as being in dense environments, in fact, most are 
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Figure 7. The kinematic misalignment angle between the ionised gas and the stars for fast-rotator galaxies, plotted against the local luminosity surface density 
within cylinders of depth 600 km s^^ in the redshift direction, with an angular size adapted to include the 3rd nearest neighbor Virgo galaxies are plotted with 
red stars, and field/group galaxies with solid blue circles. The dashed line is a guide to the eye, at the suggested ciitical density of 10^^'^ Lq Mpc^^. The error 
on each kinematic misalignment angle measurement is Ril5°. 



different. This is likely to be because the brightest and most mas- 
sive galaxies in dense environments are usually slow rotators (Paper 
VII). It seems that for fast-rotating galaxies both mass and environ- 
ment can be important independent parameters in determining the 
origin of the molecular and ionised gas. 

4 Discussion 

The results presented in the previous section shed light on the 
origin of the molecular and ionised gas in local ETGs, but they pose 
additional puzzles that require further discussion. 

4.1 The link between the gaseous phases 

Our results suggest that the molecular, atomic and ionised gas form 
a single structure, and have undergone a com mon evolution (as pre- 
viously found for ato mic and ionised gas bv iMoreanti et alj|200a : 
lOosterloo et al.ll2010n at least in fast-rotating galaxies. It would be 
difficult to maintain two kinematically distinct phases of the ISM 
for any length of time, given the collisional nature of gas (unless 
they were widely separated in radius, e.g. separate rings). For in- 
stance, let us consider a galaxy which has stellar mass loss ongo- 
ing, and has built an aligned gas reservoir. If this galaxy acquires 
misaligned gas from an external source, the two systems may in- 
teract. Collisions between two components with different angular 
momentum vectors should cause some of the gas to lose its angular 
momentum and fall to the centre resulting in a more compact gas 
distribution. A search for such a signature will be presented in a 
later work in this series. 

The ionised gas in these galaxies will also be coupled to the 
properties of the broader ISM (due to its dissipative nature). The 
dominant ionisation mechanism in ETGs is thought to be irradia- 
tion by old stellar popul a tions (e.g. post-a symptotic giant branch 
stars; iBinette et ai]|l994 ISarzi et al.ll201(]|) . which are uniformly 
distributed throughout the galaxy, hence ionised gas should be pro- 



duced at the interface of any part of the cold ISM, and the ionised 
gas would have the same angular momentum as the material it is 
produced from. 

The difference in angular resolution (and field of view) be- 
tween our CO and HI observations (and to a lesser extent our 
ionised gas observations) could have some effect on these conclu- 
sions. For instance warps are often seen in outer HI disks aroun d 
galaxies of all types (e.g. lOarcia-Ruiz. Sancisi & KuiikenI 120021) . 
However the agreement between the CO (and ionised gas) in the 
inner parts, and the HI at larger radii would be hard to manufacture 
by chance. Although challenging, higher resolution HI observations 
of the inner parts of these galaxies could be used to remove this un- 
certainty, and further establish the link between these phases. 

4.2 The importance of external accretion 

4.2.1 Slow-rotators 

In slow rotators there is clear evidence that their ionised gas is pri- 
marily acquired externally. The kinematic misalignments are statis- 
tically consistent with a flat distribution, with no significant increase 
in the number of galaxies with exactly co- or counter-rotating gas. 
This suggests that accretion and mergers dominate the gas supply 
in slow-rotating galaxies. If the source of ionization in these galax- 
ies is long lasting (e.g. irradiation by old stellar populations) then 
this gas could even be left-over from the last merger that created the 
slow-rotator. 

Slow rotators are thought to be formed from progenitors that 
have undergone multiple major mergers (Paper III). Simulations 
suggest an average of 3 m ajor mergers per slow-rota ting galaxy, 
and multiple minor mergers l lKhochfar et al. I submitted. Paper VIII). 
Slow rotators are predominately massive (Paper III), and sit at the 
centre of their group and cluster environ ments (Paper VII), whic h 
is likely to enhance their merger rate (e.g. lvan Dokkum et al.ll999n . 
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Figure 8. The kinematic misalignment angle between the ionised gas and the stars for fast-rotator galaxies, plotted against the total absolute _fCs-band magnitude 
of the galaxy from 2MASS. Distances to the galaxies are taken from Paper I. Virgo galaxies are plotted with red stars, and field/group galaxies with solid 
blue circles. The dashed line is a guide to the eye, at the suggested critical magnitude of Mx=— 24 mag. The error on each kinematic misalignment angle 
measurement is Ril5°. 
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Figure 9. As Figure^ but galaxies with Mk < —24 mag are plotted with dark green stars (for Virgo galaxies) and solid dark green circles (for field/group 
galaxies). The eiTor on each kinematic misalignment angle measurement is Kil5°. 



Our results provide additional evidence that mergers and/or accre- 
tion are important in the evolution of slow rotators. 

A few questions do remain. For instance there is no statisti- 
cally significant evidence for any of these galaxies having gener- 
ated aligned ionised gas from stellar mass loss. The stars in these 
systems are predominantly old (e.g. iKuntschner et alJuOlOl : Mc- 



Dermid et al., in prep), and hence the mass loss rate from the stellar 
population will be low. Any mass that is lost from these stars into 
the ISM may be kept hot by the large velocity dispersion of these 
pressure supported systems, and may join the X-ray halo. 

Furthermore, the lack of a peak in the misalignment distribu- 
tion for CO- and counter-rotating gas suggests the timescale required 
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Table 2. Synopsis of the statistics presented in tliis paper 



Category Number 

(1) (2) 



Mol. Detection 

(3) 



Ion Detection 

(4) 



# Mapped Mol 

(5) 



# Mapped Ion 
(6) 



Misaligned Mol 

(7) 



Misaligned Ion 
(8) 



Fast-Rotators 


224 


24±3% 


67± 3% 


36 


111 


23± 7% 


36± 5% 


Slow-Rotators 


36 


6±4% 


86± 6% 


2 


22 


50±35% 


68±10% 


Field FR 


202 


22±3% 


73± 3% 


26 


91 


28± 9% 


42± 5% 


Virgo FR 


58 


21±5% 


47± 6% 


10 


20 


10± 9% 


10± 6% 


h < lO"'' FR 


209 


22±3% 


67± 3% 


35 


93 


24± 7% 


44± 5% 


h > lO^^'^ FR 


51 


16±5% 


69±10% 


1 


18 


0% 


16± 9% 


Mks > -24 FR 


191 


24±3% 


65± 3% 


29 


92 


39± 9% 


44± 7% 


Mks < -24 FR 


33 


27±8% 


80± 7% 


7 


19 


14±13% 


5± 5% 



Notes: Column 1 lists the categories of galaxies we discuss in this paper. FR denotes fast-rotators only. I3 is our local the luminosity surface d ensity estimator, 
based o n a redshift cylinder with a depth of 600 km s^^, and an angular size adapted to include the 3rd nearest neighbor, as presented in ICappellari et al.l 
<2011bl) (Paper VII; Table 2, Column 7 and available online at http://www.purl.org/atlas3d I. M^ ^ is the absolute Ks-band total magni tude. Column 2 list s the 
number of galaxies from the ATLAS''^ survey that fall in this class. Columns 3 and 4 list the molecular gas detection rates (taken from lYoung et alj201 11) and 
ionised gas detection rates (above an integrated equivalent width of 0.02 A in either H/3 or [OIII]) for these subsamples. Columns 5 and 6 state the number of 
galaxies in this category mapped in molecular gas, and those with measurable kinematic PA's in their ionised gas maps. Columns 7 and 8 state the percentage 
of mapped galaxies which have misaligned molecular, and ionised gas respectively ('I' > 30°). 



for the ionised gas in these systems to relax into the orbital plane 
is long (especially as the ionised gas is likely to be long lived). 
This may be because these systems are in general very round (Pa- 
per III), and hence the gravitational torques will be lower than in 
disky systems. Slow rotators can also be triaxial (Paper II & Paper 
III), allowing gas in polar configurations to be stable. 

4.2.2 Fast rotators 

The molecular and ionised gas in ~36% of the galaxies in our sam- 
ple are kinematically misaligned from the stellar body, and hence 
are likely to have an external origin. There are ways in which the 
angular momentum of gas in the central regions of galaxies can 
be altered (e.g. bars whi ch produce resonant orbits and branching; 
see IWada & Habelll995h . However the kinematic alignment with 
the atomic gas at larger scales suggests this is not a dominant ef- 
fect. In addition, 36% is likely a lower limit to the importance of 
externally acquired material. Assuming that gas from mergers and 
accretion enters galaxies isotropically then one would expect on av- 
erage ~17% of galaxies to exhibit aligned kinematics (within 30°) 
by chance (30°/180° ~ 0.17), despite the gas having an external 
origin. 

The assumption that gas from mergers and accretion enters 
galaxies isotropically is a conservative one. In fact one can imag- 
ine scenarios where the gas accretion is not isotropic, but in fact 
preferential results in aligned gas. One such example would be a 
major merger that imparts a large amount of angular momentum 
to the stars in the remnant, forcing them to rotate in the plane of 
the merger. Gas that is thrown off in the merger may fall back into 
the remnant over many gigayears, and would likely fall in preferen- 
tially aligned. Some simulations have predicted that minor mergers 
may preferentially occur with the angular momentum vector of the 
accreted system aligned alon g the major axes of the dark matter 
halo (e.g. lOeason et al.ll201 ih . If such processes are ongoing, then 
we will underestimate the importance of external accretion in this 
work. How isotropic the accretion of gas from major mergers and 
cold-mode accretion onto ETGs truly is could be tested further in 
cosmological hydrodynamic simulations. 

If a large amount of external gas is present within a galaxy, the 
gas depletion timescale through star-formation will become larger 
than the relaxation time (the time it takes for gas structures to settle 
into the equatorial plane). In this case one would expect to find 50% 
of the galaxies to have co-rotating gas, and 50% to show counter- 
rotation. The lack of a large peak at a kinematic misalignment of 



180° however suggests that this criterion is seldom fulfilled, or that 
relaxation into the equatorial plane is not a symmetric process. 

IClemens et alj ( 120 lOh compared the dust production rate in 
the envelopes of evolved AGB stars with a constraint on the to- 
tal dust mass in Virgo cluster ETGs (from Herschel observations). 
They find that passively evolving ETGs (defined by having Spitzer 
mid-infrared spectra consistent with passively evolving stellar pop- 
ulations) are not detected by Herschel, and thus the dust destruc- 
tion timescale in passive ETGs must be short. They conclude that 
passive ETG's with dust are thus likely to have acquired it from 
an external source. As dust is thought to be important to shield 
nascent molecular clouds, and catalyze the conversion of HI into 
H2 molecules (e.g. lPulev & Williamsll993h . perhaps a lack of dust 
could prevent some of our sample galaxies from regenerating their 
molecular reservoirs from internal sources. 

External accretion/mergers appear to only be important in the 
field, where at least 42% of galaxies have a significant kinematic 
misalignment between the molecular or ionised gas and the stars. 
The molecular and ionised gas are nearly always consistent with a 
purely internal origin in three (partially overlapping) populations; 
fast-rotating Virgo cluster galaxies, the densest group environments 
we sample, and in the highest luminosity fast-rotators. 

As shown in Section l34l there is clear evidence that at least 6 
of the field ETGs which are kinematically aligned have gas which 
is in fact externally accreted. Including this evidence increases the 
percentage of field galaxies which have externa lly accreted gas to 
at least 46%. iKnapp. Turner & Cunniffd l ll985h have shown that 
the HI masses of early-type galaxies are unrelated to the galaxy 
luminosity, unlike in spirals. This led the m to predict that exter - 
nal sources provide all the HI in field ETGs. lMorganti et al.l ( 120061) 
andlOosterloo et al. (2010) have also shown that in almost all cases 
where HI is present in ETGs, there are some signs that it has 
been externally accreted. As it appears that the gaseous phases are 
closely linked, this suggests that the molecular and ionised gas in 
Hl-detected ETGs frequently has an external origin, even if it is 
kinematically aligned. Hence the true importance of external accre- 
tion is likely underestimated in field environments. 

If this is so, then we are forced to conclude that mergers 
and accretion must not proceed isotropically, or relaxation of mis- 
aligned gas happens over a short timescale, and preferentially in a 
prograde direction. It is also conceivable that retrograde gas sur- 
vives for a shorter length of time (due to interactions with the hot 
halo of the galaxy, or similar). Due to the different angular resolu- 
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tions (and hence spatial scales) probed by the HI and CO data how- 
ever, these conclusions require further investigation. To first order, 
one would expect the evidence for external HI accretion to become 
stronger if one were able to investigate the kinematics and morphol- 
ogy of the HI at higher spatial resolution. A more in depth treatment 
of these issues will be included in a future work in this series. 

4.3 The effect of group and cluster environments 

The striking result that galaxies in the Virgo cluster (and other dense 
environments) almost always have kinematically aligned gas sug- 
gests that dense environments affect the accretion and/or generation 
of ionised and molecular material. The challenge is thus to find a 
consistent theory that can explain the environmental dependence of 
the kinematic misalignments, while simultaneously accounting for 
the identical detection rates (and mass fractions) of the molecular 
gas in the field and clusters (Paper IV). 

Paper IV has shown that the CO-detected galaxies in Virgo 
are virialised, and hence have been in the cluster for several gi- 
gayears. The average gas depletion timescale for field galaxies in 
our CO-detected sample is «l-2 Gyr, assuming that the gas forms 
stars a t the rate predict ed by either the Schmidt-Kennicutt rela- 
tion of iKennicuttJj 1993) or a constant star formation efficiency as 



in lSigiel et al.l ( 120081) (See Alatalo et al., in prep). Recent studies 
have found little evidence that these ETGs deviate from these star- 
formation relations, and hence morphological quenching may not 
play a large role in the prol onging the lifetime of the the gas in ma- 
jority of these systems (see lShapiro et alJl20ld . lCrocker et alj|201 ll 
and Alatalo et al, in prep). 

If the ETGs present in clusters today were like current field 
galaxies when they fell in, then they should have exhausted the vast 
majority of their molecular gas by now (due to star formation) and 
should not be detectable anymore. This implies that the ETGs de- 
tected in Virgo today either must have been systematically more 
molecular gas-rich than is typical now when they fell in, or they 
must have regenerated their gas from an internal source recently. 
Of course, if ETGs were more molecular gas-rich in the past, then 
either all the gas must have been kinematically aligned (contrary 
to today's field ETGs) or it must have somehow re-aligned upon 
entering the cluster. Both possibilities however seem rather con- 
trived, so_ttieJnternalre£lenishment hypothesis seems more plau- 
sible. JSkillman et al.ll 19961) showed that spiral galaxies in the in- 
ner parts of the Virgo cluster have higher gas-phase metallicities 
than those in the field. They suggested that these systems have been 
evolving as 'closed-boxes', undiluted by metal poor gas infall, pro- 
viding some support for this hypothesis. 

It is of course possible that the molecular (and ionised) gas in 
ETGs in Virgo all comes from cooled stellar mass loss. This would 
however require that ETGs which fall into the cluster exhaust their 
pre-existing gas rapidly, so that any kinematically misaligned gas is 
destroyed, and then that the right percentage of the ETG population 
regenerate a molecular gas reservoir. 

Figure |7] suggests that high local density environments also 
result in aligned gas. We speculate that preprocessing in groups, 
before galaxies enter clusters, could perhaps be important in ex- 
plaining the kinematic alignment of the gas. As galaxies fall into 
groups, or transit through the outskirts of clusters, merger rates 
are enhanced (e.g. van Dokkum et al. 1993). Starbursts triggered 
by these violent interactions can quickly deplete the cold molec- 
ular gas, which would rapidly remove any evidence of misaligned 
material. The galaxies in Hl-poor groups and clusters would then be 
unable to re-accrete external gas as they settle into the group/cluster 



potential and become virialized. Stellar mass loss can then regener- 
ate an aligned gas reservoir over time. 

The principal difficulty of this explanation is that it requires the 
gas generation efficiency to be higher in clusters than in field galax- 
ies. Cluster galaxies have the same molecular gas detection rate and 
mass fractions as field galaxies, but all this molecular gas must be 
generated internally. About 1 1% of galaxies in the field have a sub- 
stantial molecular gas reservoir kinematically aligned with the stel- 
lar component (i.e. half of the detections), while ~20% of galaxies 
in Virgo have kinematically aligned gas (i.e. almost all of the de- 
tections). The ATLAS^'^ ionised gas detection rate is ~73% in the 
field, but only ~48% in the cluster. If this ionised gas is cooling to 
become molecular, then it must do so with a much greater efficiency 
in the cluster environment. 

Some studies of the formation of molecular clouds (e.g. 
iHeitsch & Hartmannl2008l) have suggested that they are most likely 
to form at the intersection of gas flows or shocks, such as those 
found in the spiral arms and bars of disc galaxies. Conceivably, 
if such conditions are found more often in cluster galaxies, this 
could result in a higher molecular cloud formation efficiency. The 
morphological classification of the ATLAS'"^ galaxies presented in 
Paper II reports a slightly enhanced barred fraction in the Virgo 
cluster (36±4% versus 27±2% in the field). Almost all of the CO- 
detected galaxies in Virgo and dense groups that are aligned have 
their molecular gas distributed in rings, spirals, or within bars (see 
Alatalo et al, in prep). These structures associated with galaxy reso- 
nances can perhaps provide an environment with suitable gas flows 
or shocks, that could concentrate any diffuse gas present (from stel- 
lar mass loss) and allow it to cool and form molecular clouds. The 
small difference in barred fraction within the cluster is unlikely to 
be enough to explain the apparent enhanced molecular gas forma- 
tion efficiency however. 

^Paper VII revisited the morphology-density relation JDresslen 

Il980l) for the galaxies in our sample using the kinematic classifi- 
cation presented in Paper III. It showed that the fraction of spirals 
decreases, and the fraction of fast-rotating ETGs increases as you 
approach a group or cluster core. This can be interpreted as circum- 
stantial evidence that spiral galaxies can be transformed into fast- 
rotating ETGs as they settle into a group or cluster. Therefore, it is 
possible that some of the ETGs we detect in the group and clus- 
ter environments actually fell in as spirals, which were then sub- 
sequently transformed into early-types. We would then be detect- 
ing the remnant of their star-forming gas, that would naturally be 
co-rotating (counter-rotating molecular gas is very rare in spirals). 
This explanation would however require accepting that the molec- 
ular gas mass fraction is coincidentally similar inside and outside 
the cluster. The CO-detected ETGs in Virgo are also statistically 
indistinguishable from the other ETGs in Virgo in their brightness 
at A'-band, specific angular momentum (Ar, from Paper III) , dy- 
namical mass or stellar velocity dispersion (both from Cappellari 
et al, in prep), alpha-element enhancement and stellar metallicities 
(at least in single stellar population analyses, McDermid et al., in 
prep). Nevertheless, we feel that this possibility is worthy of further 
study, as it may provide direct observational evidence of galaxies 
undergoing morphological change due to cluster environments. 

All of the mechanisms discussed above fail to satisfactorily ex- 
plain the identical detection rate and molecular gas fractions found 
inside and outside of Virgo. The low number of detected ETGs 
means that the detection rates and molecular gas fractions in the 
cluster may be identical by chance, and of course other clusters may 
be very different. Clearly observations of other clusters are required 
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to resolve this ambiguity, and pinpoint which, if any, of the above 
mechanisms are dominant. 

4.4 The effect of galaxy mass 

Figure [8] shows that the kinematic misalignments of the molecular 
and ionised gas in fast-rotating ETGs also seem to be effected by 
the mass of the host galaxy. Furthermore, this effect appears to be 
independent of the environmental effect discussed above (see Fig- 
ure |9j. The processes resulting in these kinematic alignments must 
thus act on single galaxy scales. 

Massive galaxies are better able to contain material that 
would normally be lost during energetic events (such as super- 
novae/galactic winds), perhaps resulting in a larger pool of aligned 
material that can later cool. This process, however, can only create 
aligned gas, and thus leaves us unable to explain the lack of exter- 
nally accreted gas in these massive fast-rotating galaxies. 

AGN feedback is one of the possible mechanisms which 
could explain this lack of kin ematically misaligned gas. 
ISpringel. Matteo & Hernguistl ((2003) have shown in simulations 
that galaxy mergers and accretion events that cause nuclear gas 
inflow can trigger AGN feedback, quenching star formation and 
gas accretion on a short timescale, particularly in massive galax- 
ies. AGN feedback triggered by merging events could disrupt and 
destroy misaligned gas reservoirs as they form in the most massive 
galaxies. Gas from stellar mass loss can avoid such a fate, as it will 
not flow into the centre in great quantities due to angular momen- 
tum conservation, and hence will not fuel an AGN. 

Massive ETGs are also known to host hot. X-ray emitting gas 
halos JForman. Jones & Tuckenll985h . These hot halos can heat in- 
falling cold gas clouds, preventing the gas from reaching the galaxy 
centers. In a similar way, as mentioned above, virial shocks can 
create a halo mass threshold, above which incoming gas will al- 
ways shock heat to the virial temperature, joining the hot halo 
jBimboim & Dekell l2003h . This material is virialised in the hot 
halo, and likely forced into an aligned configuration by interaction 
with the rest of the halo gas. It may then cool in an identical way to 
stellar mass loss, producing a kinematically aligned gas reservoir. 

Both of the halo effects described above are reported to oc- 
cur primarily in galaxies of mass > 10^^ Mq, a good match to 
our apparent threshold (Ma' = —24 mag, or 8x 10^" I/q.a' ). Once 
again, however, the fact that our molecular gas detection rate is in- 
dependent of mass makes explaining these results difficult. Massive 
galaxies may well retain a large reservoir of kinematically aligned 
hot gas (from internal or external sources), but in order to achieve a 
uniform detection rate as a function of mass the cooling rate from 
the halos would have to be enhanced by just the right amount to bal- 
ance the reduction due to the lack of accretion events. It is possible 
that the star formation efficiency in these systems could be lower, 
perhaps due to morphological quenching, allowing galaxies that do 
have gas to be detectable for longer. 

Slow-rotating galaxies in general have a similar mass to this 
population of aligned fast-rotators (>3xlO^'' Mq; Paper III), but 
have gas that is mainly from external sources. This suggests, as dis- 
cussed in Papers III and VIII, that slow-rotators have evolved dif- 
ferently to fast-rotators, with a greater number of external events, 
such as mergers, shaping their properties. 

5 Conclusions 

In this work we have demonstrated that a large proportion of 
the gas found in local early-type galaxies is likely supplied by ex- 
ternal processes. Gas kinematically misaligned with respect to the 
stars is common, indicating an external origin for the gas in > 36% 



of all fast-rotating ETGs. We have also shown that the ionised, 
atomic and molecular gas in local ETGs are linked, always hav- 
ing similar kinematics and thus presumably sharing a common ori- 
gin. In the field, ~42% of galaxies have kinematically misaligned 
ionised gas, confirming that mergers and accretion play a vital role 
in supplying gas to field ETGs. Slow rotators are not generally 
detected in molecular gas, but show a flat distribution of ionised 
gas kinematic misalignments, suggesting they obtain ionised gas 
primarily through mergers and/or accretion. In many cases, fast- 
rotating galaxies with kinematically aligned molecular and ionised 
gas have HI distributions that again suggest an external origin for 
the gas. This suggests that external gas sources dominate in the field 
(as has been suggested by previous authors). 

In the Virgo cluster, the molecular and ionised material in 
fast-rotators is nearly always kinematically aligned with the bulk 
of the stars, pointing to gas supplied by purely internal processes. 
Fast-rotators in dense groups also appear to always have aligned 
gas kinematics, suggesting that the local environment is important 
in understanding the environmental dependance of the gas origins. 
Given the results of Paper IV, indicating similar molecular gas de- 
tection rates and mass fractions in clusters and the field, this sets 
constraints on potential mechanisms for generating and influencing 
molecular and ionised gas in cluster environments. 

Molecular and ionised gas in the most massive fast-rotating 
galaxies also appear to always be kinematically aligned with the 
stars, independent of environment, suggesting that the alignment 
can be caused by galaxy scale processes which reduce the proba- 
bility that cold, kinematically misaligned gas can be accreted onto 
the galaxy (e.g. AGN feedback, the ability to host a hot X-ray gas 
halo, or a halo mass threshold). This also supports the picture (dis- 
cussed in Papers III and VIII) that fast-rotating galaxies (which are 
always aligned at high mass) have a different formation mechanism 
to slow-rotators (who at similarly high masses mainly have exter- 
nally supplied ionised gas). 

We tentatively suggest that preprocessing in groups may help 
to explain the environmental dichotomy in the origin of the gas. 
Merger-induced starbursts in groups and cluster outskirts could con- 
sume any kinematically misaligned molecular gas that is present, 
and leftover ionised material would be ram pressure-stripped. Once 
galaxies settle into a cluster or Hl-poor group, external gas accretion 
and mergers are suppressed, allowing stellar mass loss to regener- 
ate a kinematically aligned gas reservoir. It is possible that features 
such as bars and rings could funnel dust and gas lost by stars to 
the centre of the galaxy, or collect them together, and could explain 
the greater efficiency with which galaxies in dense environments 
must recreate their dense gas reservoirs. Alternatively, we could be 
detecting the remnant gas left over from the morphological trans- 
formation of spiral galaxies into ETGs as they enter the cluster and 
group environments. Both of these possibilities, however, fail to ex- 
plain why the detection rate of molecular gas (and the molecular gas 
mass fractions) are similar inside and outside of the Virgo cluster. 

More work to pin down the mechanism(s) creating and/or 
aligning the gas is clearly required to unambiguously determine 
which, if any, of these effects is dominant. 

It would also be beneficial to extend this sort of analysis to 
other clusters and groups, to see if the results reported hold true 
in yet denser environments. The Fornax cluster in the southern 
hemisphere and the Coma cluster in the north are obvious nearby 
targets, which should become accessible once new facilities such 
as the Large Millimeter Telescope (LMT) and the Atacama Large 
Millimeter/sub-millimeter Array (ALMA) come online. 
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Table 3. Kinematic misalignment between the ionised gas and the stars for all the detected, fast-rotating, ATLAS'^'^ early-type galaxies. 



Name 


0ion 


A</>ion 


^ ion — star 


'^^ion — star 




(deg) 


(deg) 


(deg) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


IC()676 


348.5 


26.9 


30.5 


29.6 


IC0719 


232.5 


2.5 


173.5 


26.9 


IC1024 


31.5 


8.2 


2.0 


13.6 


NGC0680 


26.0 


2.5 


26.5 


4.5 


NGC1266 


350.5 


3.5 


56.0 


7.8 


NGC2577 


88.0 


2.8 


16.0 


3.4 


NGC2764 


189.5 


5.2 


6.5 


8.6 


NGC2778 


50.5 


4.0 


5.0 


6.2 


NGC2824 


159.5 


7.5 


0.0 


8.0 


NGC2852 


188.0 


2.8 


32.0 


5.5 


NGC2859 


269.0 


1.8 


5.0 


3.5 


NGC2950 


97.5 


2.2 


16.5 


4.0 


NGC3182 


316.5 


6.0 


4.0 


9.2 


NGC3226 


36.0 


1.0 


8.0 


6.3 


NGC3245 


184.0 


4.8 


9.5 


5.6 


NGC3248 


302.5 


4.0 


179.5 


9.8 


NGC3412 


175.0 


10.5 


19.0 


12.1 


NGC3457 


339.5 


4.8 


5.5 


38.5 


NGC3499 


125.0 


4.2 


75.0 


12.7 


NGC3595 


99.5 


7.2 


99.0 


7.4 


NGC3599 


352.5 


3.2 


63.0 


17.6 


NGC3607 


302.5 


0.5 


1.0 


2.8 


NGC3610 


128.0 


3.8 


6.5 


3.8 


NGC3613 


96.0 


0.5 


3.5 


1.3 


NGC3619 


76.0 


0.5 


24.0 


3.0 


NGC3626 


165.0 


3.2 


174.5 


4.6 


NGC3648 


315.0 


5.2 


60.5 


6.0 


NGC3665 


209.5 


2.2 


4.0 


3.0 


NGC3694 


309.5 


8.0 


159.5 


12.2 


NGC3838 


137.0 


5.2 


1.5 


6.3 


NGC3941 


144.0 


1.2 


129.0 


3.7 


NGC4026 


343.0 


2.0 


18.5 


3.6 


NGC4036 


258.0 


1.5 


3.0 


1.8 


NGC4111 


197.0 


1.0 


47.5 


2.5 


NGC4119 


292.5 


18.6 


1.0 


19.7 


NGC4143 


155.0 


1.8 


165.5 


2.9 


NGC4179 


142.8 


3.1 


0.7 


4.4 


NGC4203 


198.0 


1.8 


3.5 


6.0 


NGC4251 


207.0 


1.2 


71.5 


3.2 


NGC4255 


113.0 


2.8 


2.0 


4.3 


NGC4281 


118.0 


0.5 


33.0 


1.6 


NGC4324 


239.0 


6.8 


1.0 


8.6 


NGC4382 


29.0 


9.7 


8.5 


10.8 


NGC4417 


224.0 


9.5 


4.5 


10.4 


NGC4429 


92.5 


2.0 


6.0 


3.2 


NGC4435 


198.6 


3.1 


6.1 


3.6 


NGC4473 


101.0 


9.7 


9.0 


10.5 


NGC4474 


95.5 


4.5 


16.5 


8.1 


NGC4494 


186.2 


3.1 


1.2 


6.8 


NGC4521 


12.5 


1.8 


23.5 


3.3 


NGC4596 


112.5 


0.5 


12.5 


4.5 


NGC4612 


310.3 


3.1 


17.7 


12.0 


NGC4621 


307.0 


10.0 


14.0 


11.2 


NGC4643 


68.0 


2.5 


20.0 


4.9 


NGC4684 


292.0 


7.5 


87.5 


9.0 


NGC4694 


167.0 


45.0 


157.5 


48.9 


NGC4697 


218.0 


0.5 


29.5 


2.1 


NGC4710 


208.5 


2.2 


1.0 


4.4 


NGC4753 


85.5 


1.5 


3.0 


2.9 


NGC5103 


357.0 


8.0 


38.5 


8.9 


NGC5173 


93.5 


6.0 


174.0 


17.8 


NGC5273 


71.5 


44.5 


119.0 


45.0 


NGC5353 


297.5 


0.5 


24.5 


1.1 


NGC5355 


62.0 


8.5 


33.0 


16.4 


NGC5379 


77.0 


5.0 


16.0 


11.2 
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A*ion-star 




(deg) 


(deg) 


(deg) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC5422 


348.5 


2.0 


14.5 


4.2 


NGC5485 


264.0 


4.2 


5.0 


8.0 


NGC5493 


137.5 


5.2 


16.0 


6.5 


NGC5582 


23.5 


2.5 


6.0 


3.7 


NGC5866 


126.5 


3.2 


0.0 


3.5 


NGC6014 


151.5 


10.2 


4.5 


13.3 


NGC6017 


132.5 


5.2 


0.0 


8.4 


NGC6149 


218.0 


7.5 


17.0 


8.6 


NGC6278 


337.0 


3.0 


31.5 


5.6 


NGC6798 


313.0 


1.2 


174.0 


6.9 


NGC7465 


109.5 


9.5 


57.0 


30.5 


NGC7710 


154.0 


21.2 


20.0 


46.6 


PGCO 16060 


161.5 


2.8 


2.5 


3.4 


PGC029321 


105.5 


27.2 


55.0 


46.0 


PGC035754 


47.0 


12.0 


39.0 


16.3 


PGC042549 


242.0 


5.0 


1.0 


7.8 


PGC056772 


214.0 


2.8 


23.0 


6.1 


PGC058114 


93.1 


3.1 


153.9 


10.5 


PGC061468 


73.0 


15.2 


32.0 


19.1 


UGC05408 


122.5 


10.2 


27.5 


19.9 


UGC06176 


214.0 


8.5 


13.5 


10.7 


UGC09519 


172.0 


4.5 


77.5 


6.2 



Notes: Columns 2-5 show the ionised gas kinematic PA, error, and the misalign- 
ments between the ionised gas and the main body of the stars for all the fast- 
rotating ATLAS'^'^ galaxies where a kineinatic PA was measurable. The stellar 
PA is taken directly from Paper II, and the data tables from this paper are avail- 
able to download from http://www.purl.org/atlas3d Kinematic misali gnments 
for the fast-rotating SAURON ETGs are tabulated in lSarzi et al]<2006l) . 

Table 4. Kinematic misalignment between the ionised gas and the stars for all the slow-rotating ATLAS'*'^ early-type galaxies with measurable kinematic 
misalignments. 



Name 


(/"ion 


Action 


^ion — star 


A'I'ion-star 




(deg) 


(deg) 


(deg) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 



2.8 8.5 12.3 

21.5 11.5 23.5 

15.0 176.0 18.0 

5.5 78.5 89.9 

12.2 58.0 18.8 

1.2 142.5 89.8 

5.7 4.5 7.5 
0.5 52.0 3.3 
7.5 12.5 13.7 
1.0 125.5 89.8 
1.5 139.5 25.8 
1.0 171.0 7.3 

3.8 7.5 19.4 
2.0 174.0 9.0 

8£ 4^5 12.4 

Notes: Columns defined as in Table [3J Kinematic misalignments for the slow- 

rotating SAURON ETGs are tabulated in Sarzi et al. {2006) . 



NGC0661 


227.5 


NGC1222 


54.5 


NGC1289 


276.0 


NGC3522 


192.0 


NGC3796 


67.5 


NGC4168 


102.5 


NGC4191 


178.0 


NGC4261 


103.0 


NGC4476 


219.0 


NGC4636 


32.5 


NGC4690 


191.5 


NGC5322 


102.0 


NGC5481 


233.5 


NGC5631 


305.0 


PGC028887 


216.5 
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